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! We present results from Chandra observations of 14 ultraluminous infrared 

galaxies (ULIRGs; log(L IR /L ) > 12) with redshifts between 0.04 and 0.16. The 
goals of the observations were to investigate any correlation between infrared 
color or luminosity and the properties of the X-ray emission and to attempt to 
determine whether these objects are powered by starbursts or active galactic 
nuclei (AGNs). The sample contains approximately the same number of high 
and low luminosity objects and "warm" and "cool" ULIRGs. All 14 galaxies 
were detected by Chandra. Our analysis shows that the X-ray emission of the 
two Seyfert 1 galaxies in our sample are dominated by AGN. The remaining 12 
sources are too faint for conventional spectral fitting to be applicable. Hardness 
ratios were used to estimate the spectral properties of these faint sources. The 
photon indices, T's, for our sample plus the CTiandra-observed sample from Ptak 
et al. (2003) peak in the range of 1.0-1.5, consistent with expectations for X- 
ray binaries in a starburst, an absorbed AGN, or hot bremsstrahlung from a 
starburst or AGN. The values of V for the objects in our sample classified as 
Seyferts (type 1 or 2) are larger than 2, while those classified as HII regions or 
LINERs tend to be less than 2. The hard X-ray to far-infrared ratios for the 
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12 weak sources are similar to those of starbursts, but we cannot rule out the 
possibility of absorbed, possibly Compton-thick, AGNs in some of these objects. 
Two of these faint sources were found to have X-ray counterparts to their double 
optical and infrared nuclei. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: nuclei - 
galaxies: starburst — X-rays: galaxies — infrared: galaxies 



1. Introduction 

Ultraluminous infrared galaxies (ULIRGs) are defined as galaxies with L IR = L 8 _ 1000Mm > 10 12 L 
(H =75 km s _1 Mpc -1 , q =0). Ground-based observations have shown that almost all of 
these galaxies are undergoing mergers (e.g. Sanders et al. 1988); these galactic mergers are 
thought to be the progenitors of some elliptical galaxies (e.g. Genzel et al. 2001; Veilleux et 
al. 2002) and may be a phase through which galaxies pass before a quasar is formed. ULIRGs 
in the local universe may be compared with submillimeter sources at z = 1-4 observed with 
the SCUBA instrument (e.g. Smail et al. 1997; Hughes et al. 1998). The mean properties 
(Ltr, M(fi2), and near-infrared colors) of the two classes are remarkably similar. Integration 
of the light from the ULIRG/SCUBA population shows that it may account for most or all of 
the sub millimeter/far- infrared background, as a result of the strong cosmological evolution 
of these sources. 

It is thus important to study the nature of ULIRGs at modest redshifts in order to 
understand their evolution and star formation at high redshifts. One fundamental ques- 
tion that needs to be addressed is whether the high luminosity of these galaxies results 
from starbursts or accretion onto supermassive black holes (SMBHs). Optical and infrared 
emission-line spectra suggest the energy source of ULIRGs is mostly from starbursts (e.g. 
Veilleux et al. 1997; Genzel et al. 1998) while the "warm" infrared colors of some ob- 
jects, especially the more luminous galaxies, suggest black hole driven activity (e.g. Surace 
& Sanders 1999). There may exist an evolutionary sequence of merger- induced starburst 
galaxies ("cool" ULIRGs), then "warm" ULIRGs, and then eventually Quasi-Stellar Objects 
(QSOs). If this sequence is valid, one would expect dominance by active galactic nuclei 
(AGNs) in "warm" ULIRGs and, indeed, they tend to have Seyfert-like optical and near- 
infrared spectra (see e.g. Veilleux et al. 1995, 1999a,b). Such an evolutionary sequence can 
also be tested with X-ray observations. 
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The nuclei of ULIRGs may be very heavily obscured. Therefore, observations in UV, 
optical, near-infrared and even the far-infrared may not penetrate through the dust to the 
nucleus. High resolution, high frequency radio observations can penetrate the high columns 
and are excellent probes of whether an AGN is present (Nagar et al. 2003). However, the 
bolo metric luminosity in the radio band is insignificant, and thus radio observations cannot 
prove that accretion onto a SMBH is the dominant energy source. The remaining option 
is to observe ULIRGs in hard X-rays. While the ratio of hard X-ray luminosity to infrared 
luminosity is small in nearby ULIRGs, it is not very much smaller than that of radio-quiet 
QSOs. Pure starburst galaxies, at low redshifts, do not exhibit unresolved hard X-ray (2.0- 
8.0 keV) nuclei. Starbursts, such as M82, have extended hard X-ray emission from both 
diffuse gas and X-ray binaries (Griffiths et al. 2000). At the typical distances of the 1-Jy 
sample (z ~ 0.1), the angular extent of the diffuse, hard X-ray emission in M82 would be 
~0'. / 05 (~ 100 pc) and that of the X-ray binaries would be ~0'/5 (~ 1 kpc). Starbursts 
in ULIRGs also have typical extents of <1 kpc (Soifer et al. 2000), so their hard X-ray 
emissions in our sample may be difficult to resolve with Chandra. The situation may be 
further complicated by the possible presence of a large column density of gas (N H > 1 o 23 ~ 25 
cm -2 ), which can strongly attenuate directly viewed X-rays from an AGN. If such is the case, 
X-rays emitted along the polar axis of a disk-like gaseous structure can be electron scattered 
into the line of sight, with the signature of an Fe Ka line of large (a few keV) equivalent 
width. Discovery of such Fe Ka lines may be the best determinant of an energetically 
dominant AGN. 

Previous studies of X-ray emission from ULIRGs have been made by Ptak et al. (2003) 
with Chandra and Franceschini et al. (2003) with XMM-Newton. The Ptak et al. (2003) 
sample is a volume limited sample in which the redshifts of the galaxies do not exceed 0.045. 
Franceschini et al. (2003) selected their sample from the 15 ULIRGs observed by Genzel 
et al. (1998) which included only the brightest nearby ULIRGs and only one ULIRG with 
redshift greater than 0.082. Our sample encompasses ULIRGs with greater redshifts (0.043 
< z < 0.163), and is selected to cover uniformly the IRAS color- luminosity plane. 

The organization of this paper is as follows: § 2 discusses the sample selection, § 3 the 
observations and data reduction, § 4 the analysis and results concerning X-ray structure 
and spectra, § 5 a discussion of some astrophysical consequences, and § 6 a summary of our 
conclusions. We will assume H =75 km s _1 Mpc -1 , q =0 throughout this paper. 
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2. Sample Selection 

The "1-Jy sample" of ULIRGs comprises a sample of IRA S galaxies with fluxes at 60 /zm 
exceeding 1 Jy, Lpr > 10 12 L , galactic latitude |b| > 30°, f(60 fim) > f(12 /im) (to avoid 
stars), IRAS color log(f(60 /mi) /f (100 /im)) > -0.3, and redshift < z < 0.28 (e.g. Kim & 
Sanders, 1998; Veilleux et al. 1999a,b; Kim et al. 2002; and Veilleux et al. 2002). As part of 
our sample, we selected 13 galaxies from the 1-Jy sample. Also observed were F17208-0014 
and F23365+3604 which satisfy all the criteria except the galactic latitude one (they have 
|b| < 30°), and F15250+3609 which meets all the criteria except, marginally, the luminosity 
one (it has L IR = 10 11 - 99 L ). 

The galaxies were selected to cover the full range in the key parameters Li R and 
f(25 /zm)/f(60 /im). Specifically, we have chosen galaxies that are approximately equally 
distributed over Li R and f(25 /zm)/f(60 /im) in the following 4 bins: log(Li R /L ) < 12.3, 
log(L IR /L ) > 12.3, f(25 /W)/f(60 /im) < 0.2 ("cool" ULIRGs), and f(25 /im)/f(60 /im) > 
0.2 ("warm" ULIRGs). The sample size of 16 galaxies is large enough to adequately sample 
the range of infrared luminosities and infrared colors that characterize the class of ULIRGs. 

Of these 16 galaxies, only 14 were scheduled to be observed with Chandra. Figure 1 
depicts the distribution of the entire 1-Jy sample in the log(f(25 /im)/f(60 /im)) versus 
log(Li R /L ) plane. Also indicated are those galaxies in the 1-Jy sample which have been 
previously observed with Chandra (Ptak et al. 2003), the 14 galaxies observed by us with 
Chandra, and those observed by Franceschini et al. (2003) with XMM-Newton. It is notable 
that previous Chandra observations have focused preferentially on ULIRGs with low infrared 
luminosities (log(L IR /L ) < 12.3), whereas our sample contains equal numbers of objects 
below and above 10 123 L . Furthermore, our sample contains approximately equal numbers 
in each of the four quadrants of Figure 1. This distribution allows us to test whether objects 
with certain infrared colors and luminosities are powered preferentially by stars or by AGNs. 
It is notable that essentially all objects in our sample are classified as ongoing or old mergers 
based on a comparison between the optical/near-infrared images and published numerical 
simulations of galaxy interactions (Veilleux et al. 2002). 

3. Observations and Data Reduction 

The 14 galaxies were observed between December 2002 and September 2003. Each 
galaxy was observed in a single exposure using the ACIS S3 CCD chip with the standard 
frame time of 3.2 seconds. Total exposure times, actual dates of observations, and some 
properties of the sources are summarized in Table 1. 
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Most of the data reduction and analysis was done using CIAO v2.3 with CALDB 2.23 
and XSPEC vll.2. Only a comparison of the radial profiles of two sources with models of the 
point spread function (PSF) was done using CIAO v3.0.2 and CALDB 2.25. The effects of 
the CIAO and calibration updates since v2.3 are negligible for CCD resolution observations of 
our sources due to the low signal-to-noise ratios. The data reduction followed the procedures 
outlined in the Science Analysis Threads for ACIS data on the CIAO webpages 2 . 

The position of each X-ray source was determined using the IDL routine CNTRD. The 
routine returns the X and Y positions of the centroid of a point source starting from user- 
provided initial guess positions. The R.A. and Dec were then determined using the ds9 
software from SAO based on the X and Y output of the CNTRD routine. 

Nuclear spectra were extracted for the two bright X-Ray sources F01572+0009 and 
Z11598-0112 using the CIAO tool PSEXTRACT, which creates a source spectrum, a back- 
ground spectrum, and associated response matrices. PSEXTRACT also bins output spectra 
to a specified minimum number of counts per bin. ACISABS was then applied to correct 
for the degradation in the low energy response of the ACIS chips as a result of deposition of 
contaminants on the pre-CCD filter or the CCDs. We have ignored channels below 0.5 keV 
(where the instrumental calibration is uncertain) and above 8.0 keV (where there are few 
counts) in modeling the spectra. The data were binned to both at least 15 counts per bin 
and at least 3 counts per bin. The spectra were then modeled using the XSPEC package 
(§ 4.3.1). The 15 counts per bin spectra were modeled using \ 2 statistics, while the 3 counts 
per bin spectra were modeled in c-stat mode using Poissonian statistics. Although the c-stat 
fitting approach was devised for unbinned spectra, c-stat in XSPEC performs better if the 
data are binned to at least 1 count per bin. This ensures that there are no bins with zero 
counts or any mis-match between the source and background spectra. Therefore, the data 
were binned to 3 counts per bin for the c-stat mode. For the other 12 sources, hardness ratios 
were calculated using the counts in a soft energy band (0.5-2.0 kev) and a hard energy band 
(2.0-8.0 kev), using equation 1 in § 4.3.2. The hardness ratios were then compared with 
power law and MEKAL models, photoelectrically absorbed by an intervening column. This 
method provides estimates of the spectral parameters (§ 4.3.2). Due to the low number of 
counts from most of our galaxies, we cannot place meaningful constraints on more complex 
models. 
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4. Analysis and Results 

Here, we compare the positions of the X-ray sources with the optical and near-infrared 
positions in § 4.1. Then we describe the X-ray structures in § 4.2. All 14 sources observed 
were detected with Chandra, but only two were bright enough for detailed spectral modeling 
to be performed. An analysis of the spectra of these two bright sources is presented in § 4.3.1, 
while the spectra of the others are discussed in § 4.3.2. 

4.1. Astrometry 

The positions of the X-ray peaks are offset from the infrared and optical peaks by 
typically < 1", which is consistent with the errors of measurement in the three wavebands. 
Table 2 details the positions and offsets of each source. 

4.2. X-Ray Structure 

Two (F10190+1322 and F12112+0305) of the three sources in our sample that have 
double near-infrared and optical nuclei were found to have double X-ray nuclei. The X- 
ray separations of these nuclei agree to within ~ 1" of that of their infrared counterparts. 
Figure 2 shows X-ray grey scales of F10190+1322 and F12112+0305 with infrared and optical 
contours. The infrared and optical positions of the midpoints between the two nuclei were 
shifted to match those of the X-ray midpoints. The magnitude of the R.A. and Dec shift 
applied for F10190+1322 was 1'.'2 and 0"6, respectively, in the infrared, and 1'.'42 and C/55, 
respectively, in the optical. The R.A. and Dec shift applied for F121 12+0305 was 0'.'08 and 
/ .'9, respectively, in the infrared, and O'.'O and Cf'.9, respectively, in the optical. There is 
thus weak evidence that the X-ray peaks of F10190+1322 are offset from the infrared and 
optical peaks by a fraction of an arcsecond. The X-ray peaks of F12 112+0305 are probably 
consistent with the locations of the infrared and optical peaks, after the small spatial shift 
has been applied. The linear separations of the two X-ray peaks for F10190+1322 and 
F12112+0305 are approximately 5.6 kpc and 3.7 kpc, respectively. 

The X-ray emissions of the two bright sources F01572+0009 and Z11598-0112 are con- 
centrated in the central regions of these galaxies. Figure 3 shows X-ray grey scale representa- 
tions of these Seyfert 1 galaxies with infrared and optical contours superposed. The infrared 
and optical images have been shifted so that the infrared and optical peaks match the posi- 
tions of the X-ray peaks. There is a suggestion of an E-W extension in F01572+0009 in the 
infrared and optical, as well as in the X-ray. The morphologies of the X-ray emission have 
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been investigated by comparing the radial profiles of the X-ray sources with the PSF models 
in the standard calibration library. The X-ray spectra show that most of the observed flux 
from both of these bright sources is concentrated in the range of 0.5-2.0 keV. Therefore, we 
compared the azimuthally averaged radial profiles of these sources in this energy range with 
PSFs evaluated at 1.0 keV. This comparison (Figure 4) shows that the soft X-ray emission 
is unresolved or, at best, marginally resolved. 

Most of the remaining 10 sources appear to be unresolved, with the exceptions of F16090- 
0139 and F17208-0014. F16090-0139 appears to be extended in the NW-SE direction. Its 
linear extent is approximately 8.2 kpc (3'.'0). F17208-0014 seems to be resolved with a linear 
diameter of approximately 5.2 kpc (6'.'2). The upper limits to the linear sizes of the rest of 
the sources fall in the range of 0.5-6.3 kpc. 

4.3. X-Ray Spectra 

4-3.1. The Bright Sources 

It is not surprising that F01572+0009 and Z11598-0112 are bright X-ray sources: they 
are the only type 1 Seyferts in our sample. Using the 15 counts per bin data, one can use 
X 2 statistics to evaluate models of the continuum emission of these sources. The spectra 
were first modeled with single power laws. Due to the high flux of the soft component and 
consequently the high signal-to-noise ratio in the soft energy bins, such single power law 
models underestimate the flux in the hard energy band. Therefore, a two component model 
was needed to describe the continuum spectra: a hard power law and a soft component 
represented by another power law or a MEKAL model were used. If most of the flux in the 
soft band is produced by starbursts, then the soft band flux could plausibly be represented 
by a MEKAL model (for a hot diffuse gas). The results of our modeling are listed in Table 3, 
and the spectrum and the double power law model of F01572+0009 are shown in Figure 5. 

The 3 counts per bin data were used to determine if there are weak emission lines 
in the spectra. The continua of the data were first modeled using two power laws. The 
spectral indices from these fits are consistent with those obtained from binning the data to 
at least 15 counts per bin and using x 2 statistics (see Table 3). This agreement indicates that 
binning the data to at least 3 counts per bin did not introduce any biases. The F01572+0009 
spectrum shows an excess above the power law continuum at around 6.0 keV (Figure 5), but 
this suggestion of an emission line(s) is not significant. The spectrum of Zl 1598-0112 has a 
possible emission line at an energy consistent with redshifted Fe Ka (Figure 6). We modeled 
the spectrum with a double power law as we had previously done with the 15 counts per bin 
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data. Then a narrow Gaussian feature was added to the continuum to represent the emission 
line. Using the c-stat statistics option in XSPEC, the best fit model suggests that the line 
is located at a rest energy of 7.0 keV with an equivalent width of l.O^o!? keV (Table 3). 
Emission at 7 keV in the source frame would require the iron to be highly ionized. The 
photon indices for the soft and hard power laws are 3.47 and 0.99, respectively. Figure 6 
shows the observed spectrum of Zl 1598-0112 together with the model components. 

The significance of the line cannot be tested using the F-test because the test is only 
valid for Gaussian statistics. Therefore, simulated "fake" spectra, constructed in XSPEC, 
were used to determine the likelihood that the emission line seen in Z11598-0112 is real. A 
set of 500 fake spectra were created using the FAKEIT command in XSPEC. The task uses 
the response matrices associated with the real spectra and the best-fit continuum model 
to create artificial source and background spectra. The simulated source and background 
spectra were also binned to at least 3 counts per bin. As a sanity check, we modeled the 
simulated spectra. The distribution of the photon indices in the 500 simulated spectra were 
consistent with the distribution of the photon indices for the observation of Z11598-0112, 
modeled with 3 counts per bin. Therefore, we are confident that binning the data to an 
arbitrary small number of counts per bin did not introduce any biases. We found that only 
3 of the 500 spectra showed a flux at the energy of the line exceeding the measured line flux 
minus its error bar (a conservative measure of the line flux). Thus, the line is significant at 
above the 99% level. 

It is also important to note that Zl 1598-0112 is considered to be a Narrow-Line Seyfert 1 
(NLS1) galaxy whose H/3 line width (FWHM) is 770 km s" 1 , based on the data presented in 
Veilleux et al. (1999a). NLSls tend to have steeper soft X-ray spectra than normal Seyfert 1 
galaxies, as shown independently by ASCA observations analyzed by Leighly (1999) and 
Vaughan et al. (1999). Of the 24 NLSl's studied by Leighly (1999) and Vaughan et al. 
(1999), 79% have soft flux in excess of the power law model that fits the individual spectra 
at high energies; the excess flux dominates the spectra at energies < 1.5 keV. The nominal 
power law photon indices of NLSls over the 0.6-10 keV energy band span the range of 
1.6-2.5, larger than normal Seyfert 1 galaxies (Vaughan et al. 1999). The X-ray spectral 
properties of Z11598-0112, in particular the steep soft X-ray spectrum and the flat hard 
X-ray spectrum, are consistent with a NLS1 classification. 

4.3.2. The Faint Sources 

Twelve of the fourteen galaxies that we have observed with Chandra do not have enough 
counts for the usual spectral modeling procedure. These sources have total counts in the 0.5- 
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8.0 keV band ranging from 3 to 92. In order to determine the properties of these sources, we 
used hardness ratios to estimate model parameters from XSPEC. The hardness ratio (HR) 
is defined as 

where H is the number of counts in the hard band (0.5 - 2.0 keV) and S is the number of 
counts in the soft band (2.0-8.0 keV). The hardness ratios calculated from the data may- 
then be compared with hardness ratios generated from models (such as a power law or a 
MEKAL) to determine the model parameters which describe the observations. 

Two models were assumed - a single power law and a single temperature MEKAL. For 
a single power law, the photon index (r) was varied, while the temperature (kT) was varied 
in the MEKAL model. In both models, photoelectric absorption by cold gas was included. 
For each column density (Nh) and model parameter (r or kT) pair, XSPEC generated a 
model spectrum which was then multiplied by the effective area at each energy (obtained 
from the response matrices for the actual data) and sampled appropriately. The output was 
thus a model of the number of photons detected per second as a function of energy, which 
could be compared with the observation. These simulated data were then used to calculate 
the hardness ratio as a function of N H and T or N H and kT. One can then plot contours of 
constant hardness ratio on a diagram of N H versus T (Figure 7) or N H versus kT (Figure 8). 
In each panel of Figures 7 and 8, the middle curve represents the observed hardness ratio 
and the two other curves represent the observed hardness ratio plus and minus the error. 
We have made use of 1-a values in Tables 1 and 2 of Gehrels (1986) to estimate the errors 
in our measurements. Numerical values of model parameters are listed in Table 4, and the 
footnotes describe how the errors in the hardness ratios, T, and kT were obtained. 

The reliability of this hardness ratio method can be tested by comparing its results with 
those given by the more traditional method of fitting models to the observed spectra. For this 
comparison, we used the two bright sources F01572+0009 and Z11598-0112. The hardness 
ratio method systematically underestimated the values of the photon indices compared with 
spectral fitting of single power law models when the full energy band was considered (compare 
Tables 3 and 4). This is because, as previously mentioned in § 4.3.1, a single power law does 
not adequately describe the data for these two sources; their spectra are steeper at lower 
energies (Figures 5 & 6). Tests were made in which the two methods were compared over 
narrower spectral bands (within which a power law is a good representation of the spectra) 
and the results were found to be consistent to within the errors. 

The reliability of our hardness ratio method can be further tested by comparing our 
results for the three galaxies that were also in the samples of Ptak et al. (2003) and Frances- 
chini et al. (2003) with theirs. Our results for F12112+0305, F15250+3609, and F17208-0014 



agree with these previous XMM-Newton and Chandra observations to within the errors. 
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Discussion 



5.1. The Two Bright Sources 



Based on the radial profiles of F01572+0009 and Z11598-0112 (Figure 4), the nuclear 
soft X-ray emissions are probably unresolved (F01572+0009 may be slightly extended - see 
§ 4.2). This is consistent with the X-ray emission being dominated by the Seyfert 1 nuclei. 

The spectra of these two sources cannot be described by a single power law. At least 
two components are needed: a hard power law and a soft component represented by a power 
law or a MEKAL model. These two models describe the data equally well. We find that the 
best-fit MEKAL model for the soft component in both Seyfert 1 galaxies has kT ~250 eV. 
Ptak et al. (1999) found that similar models applied to starbursts usually have a temperature 
greater than 600 eV. Therefore, the low temperatures of the Seyfert 1 galaxies suggest that 
starburst activity may not be the dominant energy source of the soft component. The same 
conclusion was drawn by Boiler et al. (2002) for F01572+0009. 

Following the analysis done on the XMM-Newton observations of F0 1572+0009 by Boiler 
et al. (2002), we can further support our claim that F01572+0009 and Z11598-0112 are AGN 
dominated through a quantitative comparison. According to Boiler & Bertoldi (1996), the 
ratio of soft X-ray (0.1-2.4 keV) to far-infrared (40-120 /im) fluxes is F SX i/F F i R ~ 10" 2 - 5 
for an unabsorbed starburst in equilibrium and 10 _1 for an unabsorbed Seyfert 1 galaxy. 
Here we use the notation SX1 for the 0.1-2.4 keV (ROSAT) band, SX for the 0.5-2.0 keV 
( Chandra) band, HX1 for the 2-10 keV band, and HX for the 2-8 keV ( Chandra) band. The 
far- infrared fluxes in the 40-120 /im band can be estimated using Equation 1 in Helou et al. 
(1985) transcribed here: 



where i v are flux densities in Jy, and Fpm is in W m~ 2 . Using the flux densities in the IRAS 
60 and 100 /im bands, we estimate the far-infrared fluxes for F01572+0009 and Z11598- 
0112 to be 9.93 x 10~ n ergs cm~ 2 s _1 and 1.13 x 10~ 10 ergs cm -2 s -1 , respectively. The 
Chandra soft X-ray fluxes are F sx = 1-54 x 10~ 12 ergs cm" 2 s" 1 for F01572+0009 and 
5.79 x 10~ 13 ergs cm -2 s _1 for Z11598-0112. The F S x values can be scaled to F S xi values 
based on the photon indices of the soft power law continuum models: 



F Fm = 1.26 x 10~ 14 x [2.58f„(60/im) + f^lOO/xm)], 



(2) 




(3) 
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where the spectral index a = T — 1 and f v oc v a . From these estimates, we find that both 
F01572+0009 and Z11598-0112 have F SX i/F F ir ^ lO" 1 - 2 . These flux ratios are approximately 
consistent with the value (10 -14 ) found by Boiler et al. (2002). The soft X-ray to far infrared 
flux ratios thus indicate that the two Seyfert 1 galaxies in our sample are energetically 
dominated by AGNs. 

Given the large obscuration to the nuclei of ULIRGs, the X-ray flux in the soft band 
may be heavily attenuated by photoelectric absorption. Therefore, a better method of de- 
termining whether a source is starburst or AGN dominated is to compare its hard X-ray flux 
to its bolometric flux. Sanders & Mirabel (1996) suggested that, on average, the bolometric 
flux (Fboi) of ULIRGs is 1.15 times the infrared flux (Fir) over the 8-1000 /im band. The 
bolometric flux of our sources can be estimated using Equation 3 of Kim & Sanders (1998) 
and the IRAS flux densities taken directly from the IRAS Faint Source Catalog: 

Fir = 1.8 x 10~ 14 x [13.48 x f v (12^in) + 5.16 x f„(25/W) + 2.58 x f„(60/W) + ^(100/W)], 

(4) 

where f u are flux densities in Jy, and Ffir is in W m^ 2 . 

A study of 109 quasars from the Palomar-Green survey by Sanders et al. (1989) suggests 
that the hard X-ray (HX1, 2-10 keV) to bolometric luminosity ratio > 10 -4 for these quasars. 
Equations "0" and 1 of Franceschini et al. (2003) suggest Fiixi/Fboi < 10 -4 for starbursts. 
Therefore, a F H xi/F bo i value > 10" 4 implies AGN dominance. F H xi/F bo i for F01572+0009 
and Z11598-0112 are 10~ 2 - 2 and 10~ 31 , respectively. This result further emphasizes that 
the two Seyfert 1 galaxies in our sample are AGN dominated. This is consistent with what 
Veilleux et al. (1999a,b) concluded based on the broad line region (BLR) luminosity to 
bolometric luminosity ratio. It appears that detection of an optical/near-infrared BLR in a 
ULIRG is a sufficient condition to predict AGN dominance in a ULIRG. 

5.2. The Twelve Faint Sources 

Applying the nominal photon index derived from the hardness ratio, we estimated the 
0.5-2.0 keV soft X-ray fluxes of the weak sources. Using Equations 2 and 3, we calculated 
the soft X-ray to far-infrared flux ratio (Fsxi/Ffir) for each galaxy in our sample, and 
found Fsxi/Ffir < 10 -3 ' 5 for all sources, well below the values quoted above for both an 
unabsorbed starburst and an unabsorbed Seyfert 1 galaxy. Furthermore, all of the weak 
sources have Fnxi/Fboi < 10~ 40 . The results are provided in Table 5. The large error bars 
in the flux ratios result from uncertainties in our estimation of the photon index. 
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5.3. Correlations with Infrared Color and Luminosity 

In Figure 9 we plot the values of Fsxi/Ffir and Fhxi/Ffir for our sample of galaxies 
as a function of the infrared luminosities and the log of the IRAS 25-to-60 fim flux ratio. It 
appears that both Seyfert l's have "warm" colors, high infrared luminosity, and high X-ray 
to far-infrared flux ratios, as expected. If having "warm" colors and high infrared luminosity 
is a pre-requisite for AGN dominance, then the other source (F12072-0444) in our sample 
in this quadrant of Figure 1, which has a low X-ray to far-infrared flux ratio, could be a 
Compton-thick AGN. 

Figure 10 shows the distribution of photon index assuming the Galactic column as a 
function of infrared color and luminosity. With the exception of one LINER, all the non- 
Seyfert galaxies have photon indices less than 2. There is no obvious correlation between T 
and infrared color or luminosity. It is evident that the Seyferts have larger Ts than the rest 
of the sample. 

If an intrinsic spectral shape is assumed, we can calculate Nh from the observed spectra. 
By fixing T at 1.7, we used the hardness ratio curves in Figure 7 to estimate Nh- The results 
are tabulated in column (8) of Table 4. None of the sources appears to be Compton-thick. 
This implies that the power source for the non-Seyfert 1 galaxies could simply be intrinsically 
weak AGNs. 



5.4. Comparison with Previous Work 

Following Figure 5 of Ptak et al. (2003), we have plotted the ratio of hard X-ray to 
far-infrared flux as a function of the IRAS 25-to-60 /zm flux ratio. We reproduced the Ptak 
et al. (2003) plot and added our results to their figure (Figure 11). The two type 1 Seyferts 
in our sample lie within the region occupied by other Seyferts and composites. The dotted 
line in Figure 11 represents the average F 2 _io kev/FFm of the pure starbursts. Our data 
agree with Ptak et al. (2003) in that the ratios of the hard X-ray to FIR fluxes of ULIRGs 
are usually similar to those of pure starbursts, suggesting that most ULIRGs are powered 
by starbursts. 

Figure 12 compares the photon indices of a single-power-law fit to our galaxy spectra 
with those from the Ptak et al. (2003) sample. This histogram shows that our spectra tend 
to have higher photon indices. The Ptak et al. (2003) single power law fits are very poor 
models of the data based on the statistical values they reported. Therefore, we have also 
compared our estimated photon indices with the photon indices from their two component 
(plasma + power law) fits (Figure 13). The histograms in Figure 13 peak at T ~1.0— 1.5 for 
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both samples. 

The far-infrared luminosity is a good measure of the star formation rate (SFR) in dusty 
systems like ULIRGs. Comparison of the SFR from X-ray measurements with the SFR from 
the FIR measurements will indicate if there is any energy contribution from sources other 
than the starburst (e.g. an AGN). If the galaxy is powered purely by a starburst, then one 
would expect its "SFR in X-rays" to equal its "SFR in FIR". There are many references 
in the literature that relate the 2-10 keV hard X-ray luminosity to SFR 2 _io kev- We have 
chosen to adopt the Ranalli et al. (2003) and Persic et al. (2004) relations because they 
appear to be the best-fits to pure starbursts (Hornschemeier et al. 2004, Figure 6b). The 
two SFR 2 -io kev values will give approximate upper and lower limits as a function of SFRfir 
3 . Figure 14 relates the SFR from the Ranalli et al. (2003) and Persic et al. (2004) relations 
to the SFR from the FIR luminosity. We have included the Ptak et al. (2003) sample in the 
plot. From the figure, it is evident that our two Seyfert 1 galaxies have X-ray luminosities 
in excess of that expected from a starburst (as discussed above). The only other source 
from our sample located above the line of equality is the LINER F04103-2838. The four 
galaxies Ptak et al. (2003) determined to have AGN contributions in their spectra (Mrk 231, 
Mrk 273, IRAS 05189-2524, and NGC 6240) are also located far above the line of equality. 
The rest of the ULIRGs are most likely powered by starbursts. 

5.5. Emission Processes 

In this section, we explore the possible emission processes that may be responsible for 
the detected X-ray emission. 

5.5.1. X-ray Binaries 

Persic & Rephaeli (2002) suggested that X-ray binaries dominate the 2-15 keV lumi- 
nosity in the absence of an AGN. There are two types of X-ray binaries. The high mass type 
(HMXB) produces X-ray emission from the accretion of wind material of an OB star onto its 
neutron star or black hole companion whereas the low mass type (LMXB) produces X-ray 



3 It should be noted here that both the Ranalli et al. (2003) and Persic et al. (2004) relations are calibrated 
based on the Kcnnicutt (1998) relation for FIR. Kcnnicutt (1998) defined the FIR band to be 8-1000/um, 
which is our definition of the IR band. However, since subsequent SFR relations are calibrated assuming the 
Kcnnicutt (1998) relation is for the wavelength range of 40-120 ^m, we will use our FIR values to determine 
SFRfir to be consistent with the literature. 
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emission from accretion onto a neutron star or black hole via Roche lobe overflow from a 
low mass companion (Persic & Rephaeli 2002). Assuming the lifetime of a single burst of 
star formation is approximately 10 s years, the X-ray emission from HMXBs are expected to 
dominate the hard spectra because the low mass stars have not had time to evolve away from 
the main sequence and to come into Roche lobe contact (Persic et al. 2004). As a result, 
the spectra of the galaxies should reflect the properties of HMXBs with average T ~ 1.0-1.4 
(Persic et al. 2004, and references therein). 

However, in mergers, multiple events of starbursts may occur due to recurrent tidal 
interactions (Persic et al. 2004). In this scenario, the low mass companions in binaries 
formed in the earlier starburst events have had time to evolve and, therefore, the X-ray 
emission from LMXBs may also contribute to the hard X-ray spectra. To determine the 
significance of this contribution, mass estimates from Veilleux et al. (2002) were used in 
conjunction with Equation 7 of Colbert et al. (2004) to approximate the SFR based on the 
global hard X-ray luminosity. The LMXB-subtracted SFRs are within ~3% of the SFRs 
found without subtraction of the contribution of LMXBs. Therefore, the LMXBs do not 
contributed significantly to the X-ray luminosity. 

Based on the peaks of the histograms in Figure 13 and the above discussion, the spectra 
of our sources may well be dominated by contributions from HMXBs. This implies that the 
weak sources are starburst dominated, as suggested by their values of log(F2_io kev/FFm) 
(Figures 11 and 14). If we assume that the luminosity in the hard X-ray band is solely due 
to X-ray binaries with luminosities > 10 37 ergs s -1 , then the galaxies in our sample contain 
10 3 -10 5 binaries that contribute to the X-ray emission. Assuming a universal stellar initial 
mass function and star formation rate, these values are 4-170 times the number of X-ray 
binaries in the nearby starburst galaxies M82 and NGC 253 which in turn have 4-16 times 
the number of binaries in the Milky Way (Persic & Rephaeli 2002). If ultraluminous X-ray 
sources (ULXs) contribute to the emission, then a smaller number of X-ray binaries would 
be required. 



5.5.2. Thermal Bremsstrahlung 

On the other hand, the hard X-ray emission may result from thermal bremsstrahlung 
from a hot wind driven by either a starburst or an AGN. For low abundances, the dominant 
emission process of a thermal gas is, of course, bremsstrahlung. Rupke et al. (2002) and 
Rupke et al.(2005a,b,c) have shown evidence for galactic winds in the 1-Jy sample. The 
spectrum of thermal bremsstrahlung is almost "flat" with V ~1.2 for E < kT. This value of 
T coincides with the peak of the histograms of observed photon indices in Figures 12 and 13. 
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For a given luminosity, a relationship between the size of the emitting region and electron 
density can be determined. Assuming the emitting region is spherical, the ion density equals 
that of the electrons, and a gas temperature of 10 7 — 10 8 K, the radius of the emitting region 
based on Equation 5.14b and Figure 5.3 of Rybicki & Lightman (1979) is 

R « 1.5 x 10 7 L|/"^nc ~ 3 cm, (5) 

where /is the filling factor for the hot gas (Loc R 3 /eff ) . If the luminosity of these sources in the 
0.5-8.0 keV band is completely due to thermal bremsstrahlung, assuming an electron density 
of n e = 1 cm" 3 and / = 1, the size of the emitting regions ranges from 0.21 to 0.50 kpc. 
These sizes of the X-ray emitting regions, deduced assuming dominance of bremsstrahlung 
emission, are consistent with our finding that most of our faint sources are unresolved by 
Chandra. 



5.5.3. Absorbed Active Galactic Nuclei 

As discussed in Section 5.3, if we assume T ~ 1.7 (the photon index of a typical unab- 
sorbed AGN), the column density can be estimated from the hardness ratios. For many of the 
galaxies in our sample, the total column density estimated in this way exceeds the Galactic 
value (see Table 4). This suggests that many of the galaxies in our sample may be absorbed 
AGNs with total N H about 2-16 times the Galactic value - i.e. up to N H ~ 5 x 10 22 cm" 2 . 
These column density estimates could be inaccurate if absorption is patchy or if scattering 
from an ionized medium is significant, as in the cases of NGC 1068 (Matt et al. 2004) and 
NGC 6240 (Ptak et al. 2003). Nevertheless, Compton-thick AGNs have been detected in 
several ULIRGs - e.g. IRAS 19254-7245 (Braito et al. 2003), Mrk 231 (Braito et al. 2004; 
Ptak et al. 2003). If the intervening column density is ~1000 or more times the Galactic 
value, then Compton scattering will become important. The quality of our data does not 
allow us to rule out this possibility. 

6. Summary 

We have obtained and analyzed X-ray observations of 14 ULIRGs with the Chandra 
X-Ray Observatory. Although all 14 galaxies were detected in the 0.5-8.0 keV energy range, 
only two were bright enough for traditional spectral fitting to be applicable. Spectral analysis 
of these two galaxies (F01572+0009 and Z11598-0112) with Seyfert 1 type optical spectra 
shows that their soft X-ray emissions are unresolved. There is a suggestion of an emission 
line at 7.0 keV in the rest frame of Z11598-0112. Monte Carlo simulations of the spectrum 
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showed that this line is significant at more than the 99% confidence level. The soft X-ray 
to far-infrared flux ratios and the hard X-ray to bolometric flux ratios of F01572+0009 and 
Z11598-0112 indicate that these two galaxies are energetically dominated by AGNs. This 
result is consistent with the statement by Veilleux et al. (1999a,b) that detection of an 
optical/near-infrared BLR in a ULIRG is a sufficient condition to predict AGN dominance 
in a ULIRG. 

Unfortunately, the rest (and majority) of our sample is too faint for conventional spectral 
fitting. Instead, we used hardness ratios to estimate their spectral parameters. The soft X- 
ray to far-infrared flux ratios and the hard X-ray to bolometric flux ratios suggest that these 
galaxies are not energetically dominated by AGNs. A comparison between star formation 
rates derived from hard X-ray and far-infrared luminosities seems to support the idea that 
these objects are powered by starbursts, although the large uncertainties on the empirical 
relation between hard X-ray luminosities and star formation rates prevent us from making a 
definitive statement. A histogram of the photon indices of the X-ray spectra peaks at T = 1.0- 
1.5, consistent with the spectrum expected from high mass X-ray binaries or bremsstrahlung 
from a hot, starburst- (or AGN-) driven wind. However, we cannot rule out the possibility 
that these very hard X-ray spectra and the very low hard X-ray to far-infrared flux ratios 
are produced by absorbed AGNs. 
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Table 1. Some Properties of the Sample and the Chandra Observations* 1 



Source 




z 


logic (^) 




Spectral 


In 1-Jy 


Observation 


Exposure 


Name b 










Type 


Sample? 


Date 


(to) 


F00188-0856 


0. 


.128 


12.33 


-0.85 


LINER 


Yes 


2003 Sep 4 


9.80 


F01004-2237 


0. 


.118 


12.24 


-0.54 


H II galaxy 


Yes 


2003 Aug 3 


9.40 


F01572+0009 


0. 


.163 


12.53 


-0.61 


Seyfert 1 


Yes 


2003 Aug 26 


10.60 


Z03521+0028 





.152 


12.45 


-1.06 


LINER 


Yes 


2002 Dec 25 


7.20 


F04103-2838 


0. 


118 


12.15 


-0.53 


LINER 


Yes 


2003 Apr 28 


10.00 


F10190+1322 





.077 


12.00 


-0.94 


H II galaxy 


Yes 


2003 Jan 31 


9.40 


Z11598-0112 


0. 


151 


12.43 


-0.64 


Seyfert 1 


Yes 


2003 Apr 13 


10.20 


F12072-0444 


0. 


129 


12.35 


-0.66 


Seyfert 2 


Yes 


2003 Feb 1 


9.20 


F12112+0305 





.073 


12.28 


-1.22 


LINER 


Yes 


2003 Apr 15 


10.00 


F15130-1958 





.109 


12.09 


-0.69 


Seyfert 2 


Yes 


2003 Jun 2 


9.80 


F15250+3609 d 





055 


11.99 


-0.74 


LINER 


No 


2003 Aug 27 


9.20 


F16090-0139 





.134 


12.49 


-1.14 


LINER 


Yes 


2003 Feb 10 


9.80 


F17208-0014 e 





.043 


12.39 


-1.27 


H II galaxy 


No 


2003 May 7 


8.60 


F23365+3604 d 





064 


12.10 


-0.94 


LINER 


No 


2003 Feb 3 


10.20 



a Redshift and IR luminosity are taken from Kim et al. (2002), unless otherwise noted. Spectral types are 
taken from Veilleux ct al. (1999a) and references therein, unless otherwise noted. 

b All source names should be preceded by IRAS. 

c Total good time interval after dead-time corrections. 

d Redshift, IR luminosity, and spectral type are taken from Surace et al. (2000) and references therein. 
c Rcdshift, IR luminosity, and spectral type are taken from Veilleux et al. (1999b). 



Table 2. Astrometry of Our Sample a 



Source RA Dec RA Dec RA Dec ARA b ADcc b ARA b ADec b 

Name X-Ray X-Ray Opt. Opt. IR IR XR - Opt. XR - Opt. XR - IR XR - IR 

(h m s) (° ' ") (h m s) (° ' ") (h m s) (° ' ") (") (") (") (") 



F00188-0856 


00 


21 


26 


.54 


-08 


39 


25.9 


00 


21 


26 


.48 


-08 


39 


27.1 


00 


21 


26 


.48 


-08 


39 


27.1 


0.90 


1.2 


0.90 


1.2 


F01004-2237 


01 


02 


49. 


99 


-22 


21 


57.3 


01 


02 


49 


.92 


-22 


21 


57.0 


01 


02 


49 


,94 


-22 


21 


57.3 


1.05 


-0.3 


0.75 


0.0 


F01572+0009 


01 


59 


50 


26 


+00 


23 


40.9 


01 


59 


50 


.22 


+00 


23 


40.6 


01 


59 


50 


.23 


+00 


23 


40.5 


0.60 


0.3 


0.45 


0.4 


Z03521+0028-E 


03 


54 


42 


22 


+00 


37 


02.9 


03 


54 


42 


.23 


+00 


37 


02.4 


03 


54 


42 


25 


+00 


37 


02.0 


-0.15 


0.5 


-0.45 


0.9 


Z03521+0028-W 
















03 


54 


42 


.16 


+00 


37 


02.4 


03 


54 


42. 


15 


+00 


37 


02.0 










F04103-2838 


04 


12 


19 


43 


-28 


30 


25.0 


04 


12 


19 


.47 


-28 


30 


24.4 


04 


12 


19 


.53 


-28 


30 


24.4 


-0.60 


-0.6 


-1.50 


-0.6 


F10190+1322-E 


10 


21 


42 


.73 


+13 


06 


55.3 


10 


21 


42 


.85 


+13 


06 


55.3 


10 


21 


42 


.81 


+ 13 


06 


55.0 


-1.80 


0.0 


-1.20 


0.3 


F10190+1322-W 


10 


21 


42. 


.48 


+13 


06 


54.2 


10 


21 


42 


.55 


+13 


06 


53.1 


10 


21 


42 


56 


+13 


06 


53.3 


-1.05 


1.1 


-1.20 


0.9 


Z11598-0112 


12 


02 


26 


.77 


-01 


29 


15.4 


12 


02 


26 


.76 


-01 


29 


15.7 


12 


02 


26 


70 


-01 


29 


15.8 


0.15 


0.3 


1.05 


0.4 


F 12072-0444 


12 


09 


45 


.15 


-05 


01 


13.8 


12 


09 


45 


.12 


-05 


01 


13.9 


12 


09 


45 


12 


-05 


01 


13.9 


0.45 


0.1 


0.45 


0.1 


F12112+0305-NE 


12 


13 


46. 


06 


+02 


48 


41.0 


12 


13 


46 


.11 


+02 


48 


42.4 


12 


13 


46 


.07 


+02 


48 


42.0 


-0.75 


-1.4 


-0.15 


-1.0 


F12112+0305-SW 


12 


13 


45 


.97 


+02 


48 


39.0 


12 


13 


45 


.92 


+02 


48 


39.4 


12 


13 


45 


97 


+02 


48 


39.8 


0.75 


-0.4 


0.00 


-0.8 


F15130-1958 


15 


15 


55. 


.20 


-20 


09 


16.9 


15 


15 


55 


.16 


-20 


09 


17.0 


15 


15 


55 


16 


-20 


09 


17.2 


0.60 


0.1 


0.60 


0.3 


F15250+3609 c 


15 


26 


59 


.45 


+35 


58 


37.1 


15 


26 


59 


.48 


+35 


58 


37.7 


15 


26 


59 


41 


+35 


58 


37.3 


-0.45 


-0.6 


0.60 


-0.2 


F16090-0139 


16 


11 


40 


43 


-01 


47 


06.4 


16 


11 


40 


.42 


-01 


47 


06.5 


16 


11 


40 


.42 


-01 


47 


05.8 


0.15 


0.1 


0.15 


-0.6 


F17208-0014 c 


17 


23 


22 


.00 


-00 


16 


59.8 


17 


23 


21 


.99 


-00 


17 


00.6 


17 


23 


21 


96 


-00 


17 


00.8 


0.15 


0.8 


0.6 


1.0 


F23365+3604 c 


23 


39 


01 


30 


+36 


21 


08.4 


23 


39 


01 


.25 


+36 


21 


09.1 


23 


39 


01 


27 


+36 


21 


08.6 


0.75 


-0.7 


0.45 


-0.2 



a Positions are given in J2000 coordinates. X-ray positions are from this work. Unless otherwise noted, optical and IR positions arc taken from Kim et al. 
(2002), who state that their typical positional error (~ 2a) is estimated to be less than 0''5. 

b ARA(Ai - A 2 ) = [RA(A 1 )-RA(A 2 )]xCos[Dcc(Ai)] x [15"/ls]. ADcc(Ai - A 2 ) = Dcc(Ai)-Dcc(A 2 ). 

c These sources are not part of the IRAS 1-Jy sample studied by Kim ct al. (2002). Optical positions are taken from the USNO A2.0 Catalog. IR positions are 
taken from the 2MASS Point Source Catalog (PSC). Both catalogs are accessible through the VizieR Service at http://vizier.u-strasbg.fr/viz-bin/VizieR. The 
expected astrometric accuracy for stars in the USNO A2.0 Catalog is 0'.'25. The astromctric accuracy of the 2MASS PSC is < Of.' 2. 
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Table 3. Spectral Models for F01572+0009 & Z11598-0112 a 

Source Model N H (Galactic) b kT Z° T d K a X 2 /dof 

Name (cm -2 ) (kcV) (Zq) 



Spectra Binned to 15 Counts per Bin 



F01572+0009 


PL 


2.6 X 10 20 




r o/i+ 01 
— r l- 2 - 4 _0.l 


K PLl 


, + 0.2 
— ' "^-O^ 


X 


10 - 


4 


175/132 




PL+PL 


2.6 X 10 20 




r 1 =2.6l° 1 1 


K PLl 


- 6 ' 8 -0.4 


X 


10" 


4 


142/131 










r 2 =o.7±°, : 4 


Kpl 2 


-3 4+ 3 ' 4 


X 


10" 


5 






PL+MBKAL 


2.6 X 10 20 


3+ - 1 


l.o r 1 =2.i±° \ 


KpL 1 


-6 2+ 03 
— °" 4 -0.4 


X 


10~ 


4 


133/130 



K M =3.1±°-* X 10- 4 



Z11598-0112 


PL 


2.3 X 


10 20 




_ r ,_3 4+U.l 




-2 4+ - 1 


X 


io~ 4 


88/58 




PL+PL 


2.3 X 


10 20 




_ r , -3 6+01 


K PLl 


_ 9 o + O.l 
— ' ! '^-0.1 


X 


io- 4 


58/57 










3+ - 1 


r 2 =o.2+° 3 2 


Kp L2 : 


-7 R+ 10 ' 4 


X 


1<)- 6 






PL+MBKAL 


2.3 X 


10 20 


1.0 ri=2.5+°J 


K PLl 


_! 7 +0.2 
- 1 - ' -0.2 


X 


IO" 4 


60/56 



K M =2.7+.g-j x IO" 4 



Source 
Name 


Model 


N H (Galactic) "Line Energy* EW r d 
(cm -2 ) (kcV) (kcV) 












c-stat B /PHA bins 


Spectra Binned to 3 Counts per Bin 


F01572+0009 


PL+PL 


2.6 X 10 20 V 1 =2.6+° \ 

r 2 =o.5_ 3 


K PLl 
Kpl 2 


— '• 1 -1.6 
-3 1+ 1S 


X 
X 


1(T 
10" 


-4 
-5 


260/252 


Z11598-0112 


PL+PL+ZGAUSS 


2.3 X 10 20 7.0+° \ l.O+J' 2 r 1= 3.5+° l 

r 2 =i.o±g-;« 


K PLl 
K PL 2 

Kli no 


_, o + O.l 
-1 5+ 1 - 5 

=3.3_ 2 ; 4 


X 
X 
X 


10" 
10" 
10~ 


-4 

•s 

6 


171/129 



a All errors are 90% confidence for each parameter. The models applied to the data arc a single power law, double power law, 
and MEKAL plus a power law. All models assume absorption by the Galactic column. Unless otherwise mentioned, all models 
were applied to data in the energy range of 0.5-8.0 keV. 

k Parameter fixed to the Galactic column, obtained from the CI AO observing toolkit accessible through 
http: / / asc.harvard.edu/toolkit/colden.jsp. 

c Mctallicity fixed at solar. 

d The photon index for the first power law (Ti) in the double power law model was determined using data in the soft band 
(0.5—2.0 kcVJ. This value was kept fixed when fitting the double power law model to the data over the whole band (0.5—8.0 
kcV). 

e The normalization of the power law models is Kpx, — photons cm 2 s 1 kcV 1 at 1 keV. The normalization of the MEKAL 
thermal plasma model is K M — - — *° — ^ J n c ii H dV where is the angular size distance (cm), and n c and n H arc the 

electron and hydrogen densities (cm -3 ), respectively. The normalization of the ZGAUSS model is Ki ine — total photons cm -2 
s~ 1 in the line. 

^This is the energy of the line in the rest frame of the galaxy. 
g Cash-statistics option in the XSPEC fitting package. 
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Table 4. Spectral Models Derived from Hardness Ratios 



Source 


N H (Galactic) a 


Total 


Hard 


Soft 


Hardness 


r d 


N H (r=i.7) c 


kT f 


Name 


(10 20 cm" 2 ) 


Counts 


Counts (H) b 


Counts (S) b 


Ratio (HR) C 




(10 20 cm- 2 ) 


(keV) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


F00188-0856 


3.21 


16.0 


6-0±i;S 


10.0+1? 


-0.25+°- 3 2 6 5 


x 1 +0.5 
1 - 1 -0.6 


36111 


79.9+^ 


F01004-2237 


1.58 


20.0 




14.011? 


-0 40+ ' 32 


1 4+0-6 


16+iL 


15-8l?2 


F01572+0009 


2.56 


4386.0 


674.01 2 ® 


3712.0+gi 


-0 69+ 02 
^"• Dy -0.02 


9 17+O.IO 
-0.07 




9 O+0.2 


Z03521+0028 


12.5 


3.0 


l-Oi 2 o'J 


2-ol 2 :I 


-0.33l ;g 


1 5+ 1 - 9 


25 +564 


12.6l^ 


F04103-2838 


2.45 


30.0 


12.01" 


I8.O+4 2 


-0.20i°;?S 


1 05+°' 35 
l.UO_g 45 


43+40 
4J_24 


79.9+^ 


F10190+1322 


3.78 


16.0 


6-0±|;« 


lo.oll? 


-0.25tg-.2S 


1 ,o+0. 50 
i - 18 -0.68 


oc+65 
o,J -30 


79.9l 73 


Z11598-0112 


2.25 


1481.0 


130.0 + ii 


1351.0+37 


-o.82l° :u i 


7+0.2 
z -'-o.i 




1 s+0.2 


F12072-0444 


3.32 


16.0 


2-0±?;l 


14.011? 


-0 75+ - 43 
u -'°-0.25 






2.0+f 7 


F12112+0305 


1.75 


51.0 


14.011? 


37.0l^ 


-0 45+ - 19 


1 5+°' 4 


ill 2 ! 


70+72 

' - y -3.9 


F15130-1958 


8.60 


38.0 


7.ol2:i 


31.0 + || 


-0 63+°' 24 


9 1 r+0.75 
z - ±o -0.65 


<17 


9+5.7 


F15250+3609 


1.56 


37.0 


5.0+H 


32.0+fl 


n 70+0.25 
u - lo -0.20 


9 97+I.24 
z - z '-0.77 


<7 


9 5 +5.4 


F16090-0139 


9.25 


27.0 


lo.oli? 


17.0+1 2 


-0.4ll°; 23 


1 C7+0.53 
J-O' -0.45 


15+ 33 


7.9ir 7 


F17208-0014 


9.96 


92.0 


30.0l|J 


62. 0l?j 


-o-35±8-i? 


1 43+0-27 

^^-o^ 


23+ 18 


H-2l 6 5 9 6 


F23365+3604 


9.36 


34.0 


14.011? 


20.01^6 


-0.l8l°;f 7 


1 10+°' 35 


50l 39 7 


79.9l 68 



a Column densities were obtained from the CIAO observing toolkit accessible through 
http:/ /asc. harvard.edu/toolkit/colden.jsp. 

b The counting errors for the faint sources (< 1000 total counts) are determined assuming Poisson statistics, using Tables 1 & 
2 of Gehrels (1986) and the total number of counts in each band. The counting errors for the bright sources are simply \/N. 

c The hardness ratio is defined in Equation 1. The errors given were determined from error propagation based on the hard 
and soft band counts (columns 4 and 5). 

d The photon index for a power law model where n(E)ocE- r . T was calculated with Njj fixed at the Galactic column density. 
The errors in the photon indices were determined from the hardness ratio limits. 

e Estimatcd total column density for r=1.7, the photon index typical of an unobscured AGN. Listed are the column densities 
required to produce the observed spectra. The column densities of several sources (F01572+0009, Z11598-0112, and F12072- 
0444) cannot be modified to reduce T to 1.7, since Njj must always be greater than or equal to the Galactic value (column 
(2))- 

f The temperature in a MEKAL model, assuming solar abundances. The model has an upper limit of kT=79.9 keV, beyond 
which the MEKAL model is indistinguishable from a thermal bremsstrahlung model. The errors were determined from the 
hardness ratio limits. 



Table 5. X-Ray and Infrared Fluxes and Luminosities 



Source F FIR F sxl a F H xi b logio(ff^) logio(^) L S xi c L H xi c L 



^bol 



Name (ergs cm 2 s 4 ) (ergs cm 2 s 4 ) (ergs cm 2 s x ) (ergs s 4 ) (ergs s x ) (ergs s x ) 



F00188-0856 


1.3±0. 


1 


X 


10" 


10 


7. 


9 +6.5 
z -6.0 


X 


10" 


15 


2 


q+2.4 
°-1.8 


X 


10" 


-14 


F01004-2237 


9.7±0 


6 


X 


10" 


11 


1 


1+8.8 
-1.1 


X 


10" 


14 


2. 


+1.7 


X 


10" 


-14 


F01572+0009 


9.9±0. 


.7 


X 


10" 


11 


5 


o+O.l 
8 -0.1 


X 


10" 


12 


1 


=+0.1 
°-0.1 


X 


10" 


-12 


Z03521+0028 


6.6±0. 


.5 


X 


10" 


11 


2 


+5.2 


X 


10" 


15 


3 


7 +9.6 
' -3.2 


X 


10" 


-15 


F04103-2838 


8.1±0. 


4 


X 


10" 


11 


1 


1+7.1 


X 


10" 


14 


3 


n + 2.5 

8 -2.2 


X 


10" 


-14 


F10190+1322 


1.8±0 


1 


X 


10" 


10 


6 


4+5-7 
-4.9 


X 


10" 


15 


1 


+ 1.7 
S -1.4 


X 


10" 


■14 


Z11598-0112 


1.1±0. 


1 


X 


10" 


10 


7. 


0+0.6 


X 


10" 


12 


2. 


g+0.1 


X 


10" 


-13 


F12072-0444 


1.1±0 


1 


X 


10" 


10 


2. 


q+17 
°-2.3 


X 


10" 


14 


4. 


0+3.6 
u -3.1 


X 


10" 


-15 


F12112+0305 


4.0±0. 


.2 


X 


10" 


10 


2 


q+1.1 
°-l.l 


X 


10" 


14 


3 


5 +i.7 

°-1.7 


X 


10" 


■14 


F15130-1958 


9.1±0 


.7 


X 


io- 


n 


3 


+1.3 


X 


io- 


14 


1 


4 +0.8 


X 


10" 


-14 


F15250+3609 


3.1±0. 


1 


X 


10" 


10 


2 


9+18 
-2.2 


X 


10" 


14 


6. 


9+5.8 
z -4.9 


X 


10" 


■15 


F16090-0139 


1.8±0 


1 


X 


10" 


10 


1 


9+11 
z -1.2 


X 


10" 


14 


1 


9+1.6 
y -1.5 


X 


10" 


-14 


F17208-0014 


1.5±0.1 x 


10" 


-9 


4. 


4+1-6 
' 4 -1.5 


X 


10" 


14 


9. 


0+3.7 


X 


10" 


■14 


F23365+3604 


3.4±0.2 x 10" 


10 


1 


1+7.0 
i -l.l 


X 


10" 


14 


3 


„+2.6 
y -2.4 


X 


10" 


-14 



-4 2+ ' 4 


-4.1+0-8 

— 0.5 


9 c+2.3 
Z -°-2.1 


X 


10 41 


8, 


q+8.4 
U -6.3 


X 


10 41 


1 


,0±0. 


.7 


X 


10 46 


_o q+3.4 


-4 2+ ' 4 
^• z -0.3 


o 4 +26 


X 


10 41 


6 


+4.9 


X 


10 41 


8 


.8±0. 


7 


X 


10 45 


, 9 + 0.001 
1 ' Z -0.001 


9 9+0.05 
z - z -0.05 


q 5+O.I 


X 


10 44 


9 


0+0.7 


X 


10 43 


1 


.5±0. 


1 


X 


10 46 


-4 5+ 1 ' 1 


-4 8+ 1 ' 1 


1 0+ 2 ' 7 


X 


10 41 


1. 


q+4.9 
y -1.6 


X 


10 41 


1 


.2±0. 


1 


X 


10 46 


-3 Q+ 2 - 8 
°- a -0.4 


a 1+0.3 


3 3 +21 


X 


10 41 


1. 


i +0.7 
^-0.7 


X 


10 42 


1 


.2±0. 


2 


X 


10 46 


a 4 +0.4 
4 - 4 -0.3 


A O + 0.4 

*-°-0.3 


7 o+6.'9 
'• a -6.0 


X 


10 40 


2. 


9+2.1 
z -1.7 


X 


10 41 


4 


.4±0. 


.2 


X 


10 45 


i 9 +0.05 
-0.05 


q 1 +0.03 

^^-O.OS 


q c+0.3 


X 


10 44 


1 


q+0.1 
°-0.1 


X 


10 43 


1 


,7±0, 


1 


X 


10 46 


o y+3.3 

°- ' -0.4 


-4 8 +0 ' 4 
4 -°-0.3 


8.2tf 2 


X 


10 41 


1. 


4+1.3 
-1.1 


X 


10 41 


9 


.9±0. 


6 


X 


10 45 


-4 2 +0 ' 2 
^• z -0.2 


, q + 0.2 




X 


10 41 


3 


+ 1.9 
8 -1.9 


X 


10 41 


8 


.2±0. 


.4 


X 


10 45 


q tr+0.2 


-4 2+ ' 2 
^• z -0.2 


7 7+ 44 
'■'-40 


X 


10 41 


3 


6 +1.9 
D -1.8 


X 


10 41 


5 


.8±0. 


.4 


X 


10 45 


_4 9 + 3 - 5 
^• z -0.4 


c 1+0.4 
°- 1 -0.4 


1 q + 6.5 


X 


10 41 


3. 


o+3.5 
S -3.0 


X 


10 40 


4 


.3±0 


2 


X 


10 45 


a 9+3.8 
^• z -0.4 


A q + 0.4 


4.8+1° 


X 


10 41 


7. 


4+6.1 

4 -5.8 


X 


10 41 


1 


.4±0. 


1 


X 


10 46 


a 5 +0.2 


_c 5+O.2 


i-s±l:i 


X 


10 41 


3 


q+1.4 
-1.3 


X 


10 41 


9 


.7±0. 


.4 


X 


10 46 


*-°-0.4 


a 9+0.3 
4 - Z -0.3 


q 9+58 
3 - z -9.2 


X 


10 40 


3. 


q+2.2 


X 


10 41 


5 


.6±0 


3 


X 


10 45 



a Soft X-ray flux in the range of 0.1-2.4 keV scaled from the 0.5-2.0 keV flux using Equation 3. The 0.5-2.0 keV flux of the 12 weak sources is based on the 
power law model derived from the hardness ratio. 

b Hard X-ray flux in the range of 2.0-10.0 keV scaled from the 2.0-8.0 keV flux calculated using an equation similar to Equation 3. The 2.0-8.0 keV flux of 
the 12 weak sources is based on the power law model derived from the hardness ratio. 

c uminosity distance calculated from Equation 2 of Kim & Sanders (1998). 
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Fig. 1. — A plot of the logarithm of the IRAS flux ratio f(25 /xm)/f(60 fim) against the 
logarithm of Ljr for the 1-Jy sample. In the top panel, asterisks represent objects observed 
with Chandra by Ptak et al. (2003), stars represent objects observed with XMM-Newton by 
Franceschini et al. (2003), and filled circles represent objects in our sample, observed with 
Chandra and reported in this paper. Note that a few objects have been observed by two or 
more groups. The bottom panel shows the distribution of optical spectral types of ULIRGs 
observed by Chandra and XMM-Newton. The upper case letters represent the spectral types 
of our sample while the lower case letters represent the spectral types of the Ptak et al. and 
Franceschini et al. samples (Si, S 2 = type 1 and 2 Seyferts, H = H II galaxies, and L = 
LINERs). The four quadrants are the four parameter bins described in § 2. Open circles in 
both panels show the remaining ULIRGs in the IRAS 1-Jy sample. 
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Fig. 2. — Linear grey scale representation of the X-ray emission from the double nuclei 
sources F10190+1322 (top) and F12112+0305 (bottom). The two left panels show X-rays 
as the grey scale with infrared K' band contours. The two right panels show X-rays as the 
grey scale with optical R band contours. The X-ray images have been smoothed to match 
the resolutions of the IR/optical images. The apparent "point sources" appearing on the 
edges of the F121 12+0305 images are artifacts from the smoothing process; each of these 
bright "point sources" corresponds to only one count and is not a real X-ray source. The 
infrared and optical images are from Kim et al. (2002). The linear separations between the 
two X-ray peaks are 5.6 kpc for F10190+1322 and 3.7 kpc for F12112+0305. 
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Fig. 3. — Logarithmic grey scale representations of the X-ray emission from the Seyfert 1 
galaxies F01572+0009 (top) and Z11598-0112 (bottom). The left panels are X-ray grey scale 
with infrared K' band contours. The right panels are X-ray grey scale with optical R band 
contours. Note the difference in spatial scale between the left and right panels. The X-ray 
images have been smoothed to match the resolutions of the IR/optical images. The infrared 
and optical images are from Kim et al. (2002). 
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Fig. 4. — Comparison of model PSFs with observed radial profiles for the two bright sources. 
The x-axis is distance in arcseconds from the centroid of the emission while the y-axis is the 
surface brightness in counts per square pixel. The error bars are errors on the net counts per 
square pixel assuming Poissonian statistics. In each panel, the histogram is the model PSF 
obtained from the PSF library at 1.0 keV. The points are observed total counts with energy 
in the range of 0.5 to 2.0 keV. This diagram shows that the soft X-ray emissions from these 
two galaxies are unresolved or, at best, marginally resolved. 
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F01 572+0009 Spectrum 




channel energy (keV) 

Fig. 5. — Spectrum of F01572+0009 with at least 3 counts per bin. The thick solid line is 
the double power law model, while the dashed lines are the two power law components of the 
model. There are hints of emission lines at around 6.0 keV, but the lines are not significant. 
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Fig. 6. — Spectrum of Z11598-0112, with at least 3 counts per bin. The thick solid line is 
the model, a combination of two power laws modeling the continuum (the dashed lines) and 
a narrow Gaussian modeling the emission feature (the dotted line). 
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Fig. 7. — Each plot shows contours of constant hardness ratio (defined in Equation 1) in the 
N H versus T plane. The middle curve represents, and is labelled with, the observed hardness 
ratio, while the other two curves represent the hardness ratios \-a away from the observed 
value (see Table 4). The horizontal dashed line represents the Galactic hydrogen column 
density, while the vertical dashed line represents T = 1.7. 
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Fig. 8. — Each plot shows contours of constant hardness ratio (defined in Equation 1) in the 
Nh versus kT plane. The middle curve represents, and is labelled with, the observed hardness 
ratio, while the other two curves represent hardness ratios 1-a away from the observed value 
(see Table 4). The horizontal dashed line represents the Galactic hydrogen column density. 
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Fig. 8 cont— 
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Fig. 9. — Plots of the log of the ratio of hard X-ray (2.0-10 keV) to far-infrared flux (calcu- 
lated from Equation 2) and the log of the ratio of soft X-ray (0.1-2.4 keV) flux to far-infrared 
flux versus the log of the 25;um to 60/im flux ratio and the log of the infrared luminosity 
between 8 and 1000/zm (the last two being the axes of Figure 1). The symbols represent 
the optical spectral classifications. Si represents type 1 Seyfert galaxies, S2 type 2 Seyfert 
galaxies, L LINERs, and H H II regions. There is a clear segregation between the two bright 
Seyfert 1 ULIRGs and the rest of our sample. The Seyfert l's have high X-ray luminosity, 
high infrared luminosity, high X-ray to infrared flux ratios, and "warm" colors. The 2-10 keV 
flux was calculated from Fhx in a similar manner to the way in which Fgxi was calculated 
from Fsx (see text, Equation 3). The dotted lines represent the divisions of the infrared 
colors and luminosity bins that were used in the selection of the sample for observation with 
Chandra (§ 2 and Figure 1). 
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Fig. 10. — Plot of the photon index (T), assuming the Galactic column, as a function of 
IRAS 25-60 /im colors and infrared luminosities. With the exception of one LINER, all 
the non-Seyfert galaxies have photon indices below 2. There is no clear correlation between 
the photon index and the infrared flux ratio, or between the photon index and the infrared 
luminosity. The plot key is the same as Figure 9 and the dotted lines represent the same 
divisions of the infrared colors and luminosity bins that were used in the selection of the 
sample for observation with Chandra (§ 2 and Figure 1). 
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Fig. 11 — Plot of log(F 2 _io kcv/F F iR) vs. log(f 25/1 m/f60 M m)- The type 1 Seyferts in our sample 
are distributed near the Seyferts, while the others are located among the starbursts and 
composites. All members of the 1-Jy sample have Log(f 2 5^ m /f6o M m) ^ -0.35 (see Figure 1). 
Here we have only included the values for the Ptak et al. (2003) ULIRGs derived from their 
global spectra. The dotted line represents the average Log(F 2 _i kev/FFm) values for the 
pure starbursts. 
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Fig. 12. — A comparison of the photon indices of a single-power-law model between our 
sample (solid hashes) and the Ptak et al. (2003) sample (dashed hashes). The boxes cor- 
responding to the two Seyfert 1 nuclei in our sample (F01572+0009 and Zl 1598-01 12) are 
indicated; their photon indices were calculated from the "hardness ratio" method (Table 4) 
as for the rest of the sample. The histogram indicates that the spectra from our sample are 
softer than those from Ptak et al. (2003). Note that Ptak et al. (2003) quoted error bars only 
for models with x 2 /dof < 1.5. Therefore, the mean error for their sample was calculated 
from models for two of the ten galaxies in their sample. 
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Fig. 13. — A comparison between the power-law photon indices of our sample (solid hashes) 
and the Ptak et al. (2003) sample (dashed hashes). The boxes corresponding to the two 
Seyfert 1 nuclei in our sample (F01572+0009 and Z11598-0112) are indicated; their single 
power law photon indices were calculated from the "hardness ratio" method (Table 4) as for 
the rest of the sample. The T's for the Ptak et al. (2003) data are their 2-10 keV results 
from the plasma plus power law models. The histogram indicates that the spectra from 
our sample are softer than those from Ptak et al. (2003). Ptak et al. (2003) quoted error 
bars only for models with x 2 /dof < 1.5. Therefore, the mean error for their sample was 
calculated from models for seven of the ten galaxies in their sample. 
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Fig. 14. — A plot of SFR 2 _i kev as a function of SFR F ir. The star formation rates are 
in solar masses per year. For each galaxy, two values of SFR 2 _ 10 kcv (connected by dotted 
lines) are plotted based on the relations from Ranalli et al. (2003) and Persic et al. (2004). 
The top value is from the Persic et al. (2004) relation, while the bottom value is from the 
Ranalli et al. (2003) relation. The solid line, represents the line of equality: SFR 2 _io kcv 
= SFRfir. It should be noted here that NGC 6240 (the leftmost pair of open squares) is 
a luminous infrared galaxy (LIRG) and is used by Ptak et al. (2003) only for comparison 
purposes because of its high hard X-ray luminosity compared to starburst ULIRGs. Plot 
key: filled squares - optically classified Seyfert 1 ULIRGs from this work, filled diamonds - 
Seyfert 2 ULIRGs from this work, filled triangles H II ULIRGs from this work, filled circles 
- LINER ULIRGs from this work, open squares - Ptak et al. (2003) AGN ULIRGs and open 
circles - Ptak et al. (2003) starburst ULIRGs. 



